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Abstract 



Vorticity dynamics of the three-dimensional incompressible Euler equations are cast into a 
quaternionic representation governed by the Lagrangian evolution of the tetrad consisting of 
the growth rate and rotation rate of the vorticity. In turn, the Lagrangian evolution of this 
tetrad is governed by another that depends on the pressure Hessian. Together these form 
the basis for a direction of vorticity theorem on Lagrangian trajectories. Moreover, in this 
representation, fluid particles carry ortho-normal frames whose Lagrangian evolution in time 
are shown to be directly related to the Frenet-Serret equations for a vortex line. The frame 
dynamics suggest an elegant Lagrangian relation regarding the pressure Hessian tetrad. The 
equations for ideal MHD are similarly considered. 

1 Introductory and historical remarks 

Hamilton's determined concentration on the idea of quaternions is often depicted by mathemat- 
ical historians as an obsession. Lord Kelvin wrote that (O'Connor & Robertson 1998) 

Quaternions came from Hamilton after his really good work had been done, and 
though beautifully ingenious, (they) have been an unmixed evil to those who have 
touched them in any way. 

Having fallen in and out of fashion over the last century and a half (Tait 1890), quaternions 
currently play an important part in the theory of 4-manifolds, through which it has been shown 
that the essential physics of particles and fields is governed by geometric principles. Fluid 
turbulence is one of the great unsolved problems of modern science. While viscosity plays a 
dominant role in the late development of an incompressible turbulent flow through the Navier- 
Stokes equations, the inviscid Euler equations determine the early and intermediate dynamics. 
The Euler fluid equations are known to be essentially geometrical, so it would not be surprising 
if quaternions were helpful in understanding their solutions. 

A quaternion can be constructed from a scalar s and a 3-vector r by forming the tetrad 1 
q = [s, r] that is defined by 



x We avoid the direct nomenclature "4-vector" because of the meaning assigned to this in gauge theories. 



q = [s, r] = si — r ■ a 
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where r ■ a = ^? 



TiGi and / is the 2x2 unit matrix. {a±, cr 2 , 03} are the Pauli spin matrices 
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(1.2) 



that obey the relations aiUj = —SijI — eijk&k ■ A multiplication rule between two tetrads 
qi = [si, 7*1] and q 2 = [s 2 , r 2 ] can easily be determined from these properties 



This shows that quaternions are not commutative, although their associativity is easily demon- 
strated. They are found to be extremely useful in modern inertial navigation systems, robotics 
& graphics that are specifically designed to control or track rapidly moving objects undergoing 
three-axis rotations (Hanson 2006, Kuipers 1999). In fact, Hamilton discovered them in the 
context of an algorithm for rotating the telescope in his observatory. If Kelvin were alive today, 
he might be forced to revise his negative opinion of their importance. 

Given the evidence, it is natural to reformulate Euler vorticity dynamics in terms of quater- 
nions, particularly in tracking a fluid particle that carries its own ortho- normal co-ordinate 
system. Instead of setting Euler variables in standard function spaces, in which delicate geomet- 
ric information might be lost, the principal aim of this paper is to investigate the Lagrangian 
evolution of these variables in appropriate quaternionic form in order to preserve their inherent 
geometric properties. The language of quaternions thus provides us with an alternative and 
unique look at the problem of Euler vortex dynamics. These manipulations are not specifically 
dependent upon the nature of the domain Del 3 but for those parts of our work where the 
direction of vorticity is discussed, the local existence in time of classical solutions is necessary 
(Kato 1972). Thus we restrict D to a three-dimensional periodic domain, although other more 
general forms of D C R 3 are also valid (see Majda & Bertozzi 2001). Otherwise our manipula- 
tions should be considered to be formal, particularly since Euler data gets rough very quickly. 

Three-dimensional Euler vorticity growth is driven by the stretching vector u ■ Vtt. This 
term plays a fundamental role in determining whether or not a singularity forms in finite time. 
Major computational studies can be found in Brachetei al. (1983, 1992); Pumir & Siggia (1990); 
Kerr (1993,2005); Grauer et al. (1998), Pelz (2001) and Hou & Li (2006). The Beale-Kato- 
Majda theorem (Beale et al. 1982) has been the main cornerstone of Euler analysis : one version 
of this theorem is the precise statement that Jq \\u)\\ (p)dT must be finite to prevent singular 
behaviour on D. A BMO-version of this theorem has been proved by Kozono and Taniuchi 

(2000) . However, it has become clear that not only the magnitude but also the direction of 
vorticity is important. The papers by Constantin (1994), Constantin et al. (1996), Cordoba & 
Fefferman (2001), Deng et al. (2004,2005) and Chae (2003,2005,2006) are variations on this 
theme. References and a more global perspective on the Euler equations can be found in the 
book by Majda and Bertozzi (2001). Shnirelman (1997) has constructed very weak solutions 
which have some realistic features but whose kinetic energy monotonically decreases in time 
and which are everywhere discontinuous and unbounded. For work on Euler limits see Brenier 
(1999, 2000) and for its dynamics in the more exotic function spaces see the papers by Tadmor 

(2001) and Chae (2003,2004). 

The new results in this paper displayed in Sections 1-6 can be summarized as follows. 

A well-known variable is the scalar growth rate a = u> ■ Scj (Constantin 1994). Associated 
with this is the 3- vector rotation or swing rate x = & x Su> , where cj is the unit vorticity and 
S = \{ui,j+Uj^) is the strain matrix. Together these form a natural tetrad 2 q = [a, x]- Theorem 

2 In Gibbon (2002) q = [a, \] was denoted as £ . The change of notation to Gothic variables for tetrads has 
been introduced to avoid confusion between these and 3-vectors. 



qi © q 2 = [sis 2 - ri • r 2 , s\r 2 + s 2 r± + n x r 2 ] . 
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1 of §2 shows that the Lagrangian advection equation for the vorticity tetrad to = [0, u] can 
then be written as 

— = q®to. (1.4) 

All these quaternionic variables are Eulerian variables; i.e., point-wise functions of space and 
time, but undergoing Lagrangian evolution in time. 

The tetrad q satisfies its own Lagrangian advection equation driven by the effect of the 
pressure Hessian P = {p t ij} through the variables a p = u ■ Pu> and x p = & x Together 
these also form a natural tetrad q p = [a p , x P ]- Figure 1 shows how Su, Pu and the three 
orthonormal vectors (£>, x, u x x) are related. In addition to (1.4), Theorem 1 also contains 
the results for the Lagrangian advection of q and q p . Simply stated this is 

^ + q©q + q p = 0. (1.5) 

The result in (1.5) enables us to prove a Theorem in §2 on the direction of vorticity based on 
Lagrangian trajectories: "Provided \\Xp\\l°°(v) ^ s integrable in time up to t* > on a periodic 
domain ID, no Euler singularity is possible at t* , with the exception of the case where u becomes 
collinear with an eigenvector of P at t*" . Although different in detail, this result is in the same 
style as the direction of vorticity theorems cited above and is directly a variation of the BKM 
theorem. Ohkitani and Kishiba (1995) have observed in computations that at maximum points 
of enstrophy, u: becomes collinear with the most negative eigenvector of P. Collinearity may 
therefore be an important process in vorticity growth. The pressure Hessian P and its interplay 
with the strain matrix S has appeared in the Euler and Navier-Stokes literature in various places; 
see the references in Galanti et al. (1997), Majda and Bertozzi (2001) and Chae (2006). 

At each point in space-time a fluid particle carries its own ortho-normal co-ordinate system 
(o>, Xi & x X) '■ see Figure 1. Explicit equations for Lagrangian time derivatives of this frame are 
given in §3. The corresponding Darboux vector is the particle rotation rate. The frame-equations 
are then shown to be directly related to the Frenet-Serret relations of differential geometry that 
govern the curvature and torsion of a vortex line through the arc-length derivative of its tangent, 
principal unit normal and bi-normal. Using Ertel's Theorem, explicit differential equations for 
the curvature and torsion are then found. 

It is shown in §4 how to find Lagrangian differential equations for a p and x P - The relation 
between q p and q is given in Theorem 3 where they are shown to satisfy 

-^ = q©q P + ^- (1-6) 

is a tetrad linear in q and q p whose arbitrary scalar coefficients, in principle, are determined 
by the Poisson pressure relation. 

The vorticity vector-field u ■ V is frozen into the Euler flow. Any system with a frozen-in 
vector-field will also have an associated form of Ertel's Theorem, and a corresponding tetrad 
q = [a, x]- Thus the Lagrangian-quaternionic format displayed in this paper is more generally 
applicable, as illustrated by the equations for ideal MHD in §5. Two time-clocks and two tetrads 
q ± = [a^, x^ 1 ] appear as a result because of the two Lagrangian derivatives that naturally arise 
through the use of Elsasser variables. 

Previous attempts at formulating Euler vorticity dynamics using quaternions have met with 
only partial success. Past results have appeared in reverse order : the relations between a and 
X to be displayed in Theorem 1 were derived first by Galanti et al. (1997) (see also Gibbon 
et al. 2000), which were then shown to be expressible in a quaternionic form (Gibbon 2002). 
That story was incomplete, however, because the Lagrangian advection equation for to was 
missing, as were the ideas on particle frame dynamics, the pressure relation (1.6), and results 



22nd June 06: http://arxiv.org/abs/nlin.CD/0512034 Quaternions & particle dynamics. . . 4 



on the direction of vorticity. Roubtsov & Roulstone (1997, 2001) have also formulated semi- 
geostrophic theory in terms of quaternions. 

2 Vorticity dynamics in quaternion form 




w x x 



X 




Figure 1: A vortex line with unit tangent vorticity vector u). The normal vectors x = & x Su> and 
X p = Cjy.P£j are defined in (2.9) & (2.10). Thus the three unit vectors [<2i, %, £) x %] form an ortho-normal 
co-ordinate system. Moreover, u>, Su: and u) x x are co-planar, as are cl>, Pw and u> x \ P - 



In their vorticity form, the three-dimensional incompressible Euler equations are 

Du 

u> ■ Vm = 5a; , 



Dt 



(2.1) 



where the strain matrix is written as S = \ (y,ij + uj^) and u = cmlu is the vorticity (Majda 
and Bertozzi 2001). Equation (2.1) arises from taking the curl of the Euler equations in their 
velocity formulation 

^ = _Vp, divu = 0, (2.2) 

Dt J 

in which the Lagrangian (material) derivative is defined as 

D _ d 
~Dt~ dt 

The vorticity can be expressed as a tetrad by taking the quaternionic curl of il = [0, u] 



+ u ■ V 



(2.3) 



V © il = [— diva, curlw] 



(2.4) 



Thus there exists a natural vorticity tetrad to has the divergence-free constraint built into it 



[0, w] 



(2.5) 



The results in this paper employ Ertel's theorem, which is widely used in geophysical fluid 
dynamics in the study of potential vorticity: see Hide (1983, 2004) and Hoskins, et al. (1985). 
More generally it applies to any fluid system whose flow preserves a vector field, as the Euler 
equations preserve u> • V. For the extensive history behind this result, which seems to have 
originated with Cauchy, see Truesdell & Toupin (1960), Kuznetsov and Zakharov (1997) and 
Viudez (1999). The most general form of Ertel's Theorem says that if a? satisfies (2.1) then for 
an arbitrary differentiable vector 9 



(2.6) 
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The choice of 9 as the Euler velocity field u (Ohkitani 1993) implies that the vortex stretching 
vector uj • V« = Su is governed by 

where P = {p,ij} = {d 2 p/dxidxj} is the Hessian matrix of the pressure. Thus the combination 
of (2.1) and (2.7) gives Ohkitani's relation (Ohkitani 1993) 

D 2 u 

— + Pu, = Q. (2.8) 

To understand how the direction in which the vorticity vector stretches (compresses) in relation 
to its growth rate requires an understanding of its relationship with the matrices S and P. The 
scalar and vector variables a and x are defined by 

a = u ■ SCb , x = <*> x Su> , (2.9) 

a p = u ■ Pu , x P = & x P&> ■ (2-10) 

The left part of Figure 1, based upon Su), shows the ortho- normal co-ordinate system u>. x an d 
u x x\ the right hand part of the figure shows the same figure with S replaced by P. Thus Su 
can be decomposed into its components along the two orthogonal vectors u> and x x u 

Su) = au> + x x &- (2-11) 

From (2.1) and (2.11), the Lagrangian derivatives of \u)\ and u> are given by 

D\u>\ , Du . 

— ■ — ■ = a L) , = v x uj . (2.12) 

Dt 1 '' Dt A V ; 

The quantities (a, x) are respectively the rates of change in vorticity magnitude and direction; 
that is, one may call respectively call a and x the stretching rate 3 and the rotation or swing 
rate. These variables form natural tetrads associated with to = [0, u>] 

q = [a,x] , % = [<*p,X P ] ■ (2-13) 

The following theorem shows how Euler vorticity dynamics can be formulated using quaternions. 

Theorem 1 [Euler vorticity dynamics in terms of quaternions:] The vorticity tetrad 
to(x,t) satisfies the relation 

— = q®VO, (2.14) 

while Ohkitani's relation (2.8) becomes 

D 2 tv , , 

_ 5 - + q p ©tt) = 0. (2.15) 

The tetrad q(x,t) satisfies the compatibility relation (Riccati equation) 



Dq 

Dt 



^7 + q©q + q P = o. (2.16) 



3 a and a p are Rayleigh quotient estimates for eigenvalues of S and P respectively although they are only exact 
eigenvalues when w aligns with one of their eigenvectors. Constantin (1994) has a Biot-Savart formula for a. 
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Remark: In terms of a and x, the components of (2.16) were originally calculated by an 
indirect route in Gibbon (2002), although at that time (2.14) was not yet available. Moreover 
the present formulation simplifies the proof. 

Proof: (2.14) follows from (2.1) and (2.11) by direct calculation 

— = [0, aw + X x w] = [a , x] ® [0, w] = q © hJ . (2.17) 
Following (2.11) and Figure 1, we have 

Pu = a p uj + XpXu =>• [0, Pu)] = q p © tt> . (2.18) 
Consequently, Ohkitani's relation (2.8) implies 

= -^[0, Su] = -[0, Pu] = -q p ®to, (2.19) 
which is (2.15). Differentiating (2.14) again and using (2.15) gives the compatibility relation 

^©ttj + q©(q©tt)) + q p ©tt) = 0. (2.20) 

The result (2.16) in Theorem 1 follows because of the associativity property. ■ 

The meaning of x now becomes clear. For structures such as straight vortex tubes or flat 
sheets, u aligns with an eigenvector of S and thus x = 0, in which case a is an exact eigenvalue 
of S. The Ricatti equation for q in (2.16) reduces to a simple scalar form. However as soon as 
a tube or sheet bends, twists or tangles, % / and the full tetrad is restored. Because all our 
variables are functions of (a? , t), equations (2.14) and (2.16) govern the vorticity dynamics at all 
points and all times in the flow provided solutions remain finite. 

The BKM-theorem (Bealeei al. 1984) shows that the time integral J*q ||^||l°°(d) dr must be 
finite at a time t* to rule out singular behaviour. Variations on this theme are the direction 
of vorticity theorems expressed in the work of Constantin et al. (1996), Cordoba &; Fefferman 
(2001), Deng et al. (2004,2005) and Chae (2006). Chae's result (his Theorem 5.1) is based 
on control of the time integral of \\Su> ■ Plj||oo, which is derivable from (2.7). Here, a direct 
consequence of Theorem 1 concerns the pressure Hessian and its associated variable x P ■ 

Theorem 2 : On the domain D = [0, L]p er there exists a global solution of the Euler equations, 
u G C([0, oo];H s ) n C 1 ([0, oo]; if 8-1 ) for s > 3 if for every t* > 

||xJU°°(B)dT < oo, (2.21) 



excepting the case where ui becomes collinear with an eigenvector of P at t* . 

Remark: The theorem does not imply that blow-up occurs when collinearity does; it simply implies 
that under condition (2.21) it is the only situation when it can happen. Ohkitani (1993) and Ohkitani 
and Kishiba (1995) have noted the collinearity mentioned above; they observed in Euler computations 
that at maximum points of enstrophy, ui tends to align with the eigenvector corresponding to the most 
negative eigenvalue of P. 

Proof: Consider the relation for q in Theorem 1 in (2.16) expressed in a~x components 

Da 9 9 Dy „ . 

= X 2 - o? - a p = -2a X - X • X P (2.22) 
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where a = a(x, t) and x(t) = x( x > 0- From (2.11) we know that 

\SCj\ 2 =a 2 + X 2 (2.23) 

and so 

— — <-aiuH 9N1/9 ■ (2.24) 

Dt ' 1 {o?+x)l v ; 

Because £)|u;|/-Dt = a|w| from (2.12), the magnitude of vorticity \uj\ cannot blow-up for a < 0. Thus 
our concern is with a > 0. Inequality (2.24) becomes 

^|^<KI + |X P I- (2-25) 
Integrating along bounded Lagrangian trajectories X(t,Xo) that satisfy X t = u(X(t,Xo),t), we have 

Su(X(t,x o ),t)=Su(X(0,x ),0) + f (\a p (X(T,x ),T)\ + \x p (X(T,x ),T)\) dr. (2.26) 

Jo 

When there is no collinearity between u) and Pu), the assumption of point-wise in space integrability in 
time of x p m (2-21) also extends to a p . Thus ||S , w|| tx) is bounded if (2.21) holds and, in consequence, 
so is |co | . The Bcalc-Kato-Majda theorem then guarantees regularity of the Euler equations. However, 
there still exists the possibility that \PCj\ could blow up simultaneously as the angle between ui and 
PCj approaches zero while keeping x p finite; under these circumstances J ||Xpl|L oo (D)^ T < oo, whereas 
Jo ll a plU 00 (D)^ T — * thus blow-up would still be theoretically possible. ■ 



3 Lagrangian frame dynamics for particles and the Prenet-Serret equations 
3.1 Frame dynamics for particles 

The Lagrangian dynamics of the ortho-normal frame (d>, x, ui x x) can now be evaluated. Figure 2 
illustrates the motion of a particle from one co-ordinate point in space-time to another 




Figure 2: Vortex lines at two different times t\ and with the dotted line representing the particle (•) 
trajectory moving from [x\, ti) to (x2, t%). The orientation of the ortho-normal unit vectors (uj, x, wx£) 
is shown at each of the two space-time points. 

To find a closed expression for the Lagrangian time derivatives of the ortho-normal set (o>, Xi & x x) 
requires the derivative of X- To find this it is necessary to use the fact that the 3-vector PCj can be 
expressed in this ortho-normal frame as the linear combination 

PCj = a p Cj + cix + c 2 {Cj x x) ■ (3-1) 

where the coefficients C\ and c 2 are defined by 



ci = u> ■ (% x x P ) 



C2 = -(x- X P ) 



(3.2) 
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The 3-vector product u) x PCj yields 

X p = ci(w x x) - C2X- (3.3) 

The Lagrangian time derivative of x comes from the 3-vector part of equation (2.16) for the tetrad 
q = [a, x] m Theorem 1 



Dy Dy 

- I ± = -2ax~X P -^ = -2a X + c 2 , (3.4) 



where \— \x\- Using (3.3) and (3.4) there follows 

D X -If. »\ £>(^ X X.) - , -1 - /q r\ 

— = -c« >xx), jj- =\^ + cix X- (3.5) 

Thus, according to Euler's fluid equations, the Lagrangian time derivatives of the ortho-normal set can 
be expressed as 

DCj 

= V x Q (3.6) 

Dt v ' 

= Px(^xx) (3.7) 

f = (3-8) 
where the "Darboux angular velocity vector" T> for the ortho-normal frame is defined as 

V = x --u> with \v\ 2 = x 2 + %. (3Q) 

X X 

3.2 Frame dynamics and the Frenet-Serret equations 




Figure 3: The ortho-normal frame (u>, u: x x, X.) as the Frenet-Serret frame. 

With u> as the unit tangent vector, x as the unit bi-normal b and u> x x as the unit principal normal n, 
the matrix F can be formed 

F=(i T ,(ixx) T , f), (3.10) 

and (3.6)-(3.8) can be re-written as 

OF ( -X 

— = AF, A= X -dx" 1 I . (3.11) 



Dt 








For a space curve parameterized by arc-length s, then the Frenet-Serret equations relating dF/ds to the 
curvature n and the torsion r of a vortex line are 



dF 

— = BF where B = | -k Or). (3.12) 
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It is now possible to relate the t and s derivatives of F given in (3.11) and (3.12). At any time t the 
integral curves of the vorticity vector field define a space-curve through each point x: these space curves 
are called 'vortex lines'. The arc-length derivative d/ds is defined by 

— =w-V. (3.13) 
ds 

The evolution of the curvature k and torsion r of a vortex line may be obtained from Ertel's theorem in 
(2.6), expressed as the commutation of operators 



d_ 

ds 



D d 
Dt' ds 



0. (3.14) 



Applying this to F and using the relations (3.11) and (3.12) gives the Lax pair 

aB + ^§-^ + [B,A}=0. (3.15) 
Dt ds 

Thus Ertel's Theorem gives explicit evolution equations for the curvature n and torsion r that lie within 
the matrix B and relates them to c\ and \. 

Finally we remark that the frame dynamics along each Lagrangian trajectory may be characterized 
by a curve in the M 2 or C 1 plane. For example, one might consider the quantity 

1>(s,t) = \V\exp(iJ [ax'^MdA (3.16) 

evaluated along each Lagrangian trajectory. This complex representation of the Darboux vector's effect 
is reminiscent of the Hasimoto transformation 

ip(s,t) = K(s,t)expfi J r(s',i)ds'^ (3.17) 

used for representing the propagation of a Kelvin wave along a vortex filament in terms of its induced 
curvature and torsion (Hasimoto 1972). Such a representation is potentially useful as a diagnostic for 
characterizing frame dynamics in an experimental or computational fluid flow. Thus, because the Darboux 
vector has only two components, a representation exists that reduces the description of frame rotation 
for each fluid element to a curve in a plane. 



4 A Lagrangian advection equation for q p 

One of the hurdles in pursuing a Lagrangian approach to the Euler equations is the problem of the non- 
locality of the pressure field. Overtly, we have no Lagrangian differential equations for either a p or x P : 
the usual numerical procedure is to up-date the pressure through its Poisson equation — Ap = u^jUj^. 
How to address this issue can be illustrated by an example. Differentiating the orthogonality relation 
X ■ = and using the derivative of Cj in (2.12) gives 

U) ■ =0 => = On 4.1 

Dt Dt H v ; 

where qo lies in the plane perpendicular to a) in which \ an d X P a l so ue - Thus q = A(x, t) x+B(x, t) x p 
independently it is known that A = —2a and B = —1 from (2.16). The same method may be used for 
X p by differentiating x p ■ d> = 0, leading to 

Dx 

~£H~ = X X Xp + q Where 9 = MX + A Xp (4.2) 

where n = n(x,t) and A = X(x,t) are unknown scalars. Explicitly differentiating x P = & x P& m (4-2) 
gives 

D(Pu>) 

u> (X ■ Po>) - a pX + A x x m ' = u> (x ■ Pu>) + q ■ (4.3) 
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Using the cross product x = o> x So), this can be manipulated into 

D{Puj) 
Dt 



-a p Su:}=q, (4.4) 



which means that 

— — — = a p bu) + q x uj + eui (4.5) 
where e = e(x, t) is a third unknown scalar in addition to fi and A in (4.2). Thus the Lagrangian derivative 

Dar, 



of a p = Q> ■ Puj is 



Jap + x-Xp + e- (4.6) 

Lagrangian differential relations have now been found for x P an d <Xp, but at the price of introducing 
the triplet of unknown coefficients /i, A, and e which must adjust in a flow to take the Poisson pressure 
constraint into account - they cannot be regarded as arbitrary. 

Dimensional analysis on the various Euler variables governed by equations 4 (4.6) for a p and (4.2) for 
X p shows that [w] = T -1 , [a] = T _ \ [x] = T' 1 , whereas [a p ] = T~ 2 , [ Xp ] = T~ 2 . This means that 
[A] = T -1 , [fj] = T~ 2 and [e] = T -3 . Since the Euler equations possess no other time scale fi, A and e 
must be expressible in terms of these units or their ratios 

[i = a, X , "pjXp) such that \p] = T~ 2 , (4.7) 

A = X(lu, a, Xi a v> Xp) such that [A] = T~ 1 , (4.8) 
e = e(uj, a, X , a p , Xp ) such that [e] = T -3 . (4.9) 

Now re-define the triplet such that 

A = a + Ai, n = a p + (j,i, e = -2% ■ x P + Mi" + ^i a p + £ i (4-10) 

where the new triplet is subsumed into the tetrad (the unit tetrad is I = [1, 0]) 

«P = /iiq + A 1 q p + e 1 I. (4.11) 



Theorem 3 [Dynamics of q p ] : The pressure tetrad q p = [a p , x p ] satisfies 

©q P + qj, (4.12) 



Dc\ p 



Dt 

where the triplet of scalar variables fii, Ai, and E\ within ?${x ,t) is determined by the Poisson equation 
for the pressure 

-TrP = TrS 2 - \u? . (4.13) 

Remark: It has yet to be understood what effect the Poisson pressure constraint has on the triplet of 
scalars fix, Ai, and E\ within ^J. They are not all likely to be zero because, for example, in the simple 
case of Burgers vortex a = ao — const ; a p = — a§ ; X — X p = £i = ajj and Ai = /iti = 0. 

5 Ideal MHD 

As already indicated, these ideas can be pursued for other systems that possess vortex stretching. The 
equations of ideal incompressible MHD couple an ideal fluid to a magnetic field B 

Dii 

— = B-VB-V P , (5.1) 
DB 

— =B-Vu, (5.2) 



i S and P are subsumed into the scalars a, a p , \ an d Xp so there is no need to consider them separately. 



22nd June 06: http://arxiv.org/abs/nlin.CD/0512034 Quaternions & particle dynamics. . . 11 



together with divtt = and divB = 0. The pressure p in (5.1) is the combination p = pf + i_B 2 where 
Pf is the fluid pressure. Elsasser variables are defined by combining the u and B fields such that 

v ± =u±B. (5.3) 

The existence of two velocities v ± means that there are two material derivatives (and two time clocks) 

In terms of these, (5.1) and (5.2) can be rewritten as 

= -Vp, (5.5) 



Dt 

with the magnetic field B satisfying 

D ± B 



Dt 



B ■ Vv ± = a ± , (5.6) 



together with divw* — 0. The cr ± -strctching vectors defined in (5.6) obey an Ertel's relation already 
proved in Gibbon (2002) 

^ = ~PB. (5.7) 
Dt K ' 

The relations in (5.6) thus allow us to define 

a ± = B-(B-Vv ± ) X ± = Bx{B-Vv ± ) (5.8) 

having used Moffatt's analogy between the vectors u> and B (Moffatt 1978). Moreover, because <r ± 
defined in (5.6) lie in the plane of the unit vectors B and B x x we have the decomposition 

er* =a ± B + x ± X B. (5.9) 

Thus it is easy to prove that 

which are the equivalent of (2.12) for the Euler equations. The play the role(s) of scalar Elsasser 
stretching rates, with x ± as the rotation rates. One may also define corresponding variables based upon 
the Hessian matrix P 

a pB = B-PB, x pB = BxPB. (5.11) 
We define the tetrads q ± and q pB as follows 

W B = [0, B] , q ± = [a ± , x ± ] , q pB = [a P B , X pB ] ■ (5-12) 

Theorem 4 The magnetic field tetrad to b satisfies the two relations 



D^b 



Dt 

D^ fD ± tv, 



Dt V Dt 

The tetrads q ± satisfies the compatibility relation 

D^q± 



q^tus, (5.13) 
^q pB ®tv B = 0. (5.14) 



Dt 



q ± ©q T + q pB =0. (5.15) 
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Proof: The proof of (5.13) follows immediately from (5.6) 

D ± tv B 



Dt 



D ± B 

0, 



Dt 



= [0,tr ± ]=q ± ®to B , (5.16) 



where we have used (5.9) at the last step. The proof of (5.14) follows by combining (5.6) and (5.7) 
together with the fact that PB lies in the plane of the unit vectors B and B x x p b- Thus we have the 
decomposition 

PB = a pB B + XpB xB, => [0, PB] = q pB © B . (5.17) 

The proof of (5.15) follows as a compatibility relation between (5.13) and (5.14). ■ 

Finally MHD-Lagrangian frame dynamics, in the spirit of §3, needs to be interpreted in terms of two 
sets of ortho-normal vectors B, x ± , (B x j^) acted on by their opposite Lagrangian time derivatives. 
After some calculation we find the equivalent of (3.6)-(3.8) and (3.9) is 

D T B 

— = D^xB, (5.18) 

^(Bxx*) = V^xiBx^), (5.19) 

= V^x± ± , (5.20) 
where the pair of Elsasser Darboux vectors T>^ are defined as 

V^ = x T -^B, cf = B-lx ± x(x pB + a ± X T )}- (5.21) 

6 Summary 

The tetrad reformulation of Euler's equations in this paper appears to be completely natural, giving 
results that are remarkably simple in their expression. It also provides a new hybrid picture of ideal 
fluid dynamics in which the Lagrangian fluid parcels carry ortho-normal frames, whose rotation velocity 
depends on the local Eulerian values of the pressure and vorticity. These frames are defined by three 
ortho-normal vectors: (a) along the vorticity; (b) along its rate of change following the Lagrangian 
trajectory; and (c) along the cross product of these two unit vectors. This frame is governed by the 
Darboux vector that has components that lie only in the x-c2> plane. 

Remarkably, a picture similar to that for Euler fluids also emerges for magnetic fluids described by the 
ideal MHD equations. The MHD equations have two characteristic velocities, corresponding to the two 
Elsasser variables. Thus, MHD-Elsasscr variables summon two Lagrangian characteristics along which 
the evolutionary equations reduce to ortho-normal frame dynamics. Instead of being attached to the 
vorticity vector, both of these MHD frames are attached to the magnetic field vector (Moffatt 1978). 
The second vector in each moving frame is obtained by the rate of change of magnetic field along the 
Elsasser characteristic. The two frames are then completed by taking the cross product of the first two 
unit vectors in each frame. Again the rates of rotation of these Elsasser frames depends on local Eulerian 
properties and the respective Darboux angular velocity vectors have only two components in each frame. 

An interesting direction of future work would be to numerically monitor the tetrads q and q p to see 
how close the relations between them are adhered to. To remind the reader of the relation between them, 
we proved in Theorem 1 that q satisfies 

^ + q©q + q p -0, (6.1) 
and in Theorem 3 it was shown that if the triplet of scalars is chosen in a certain way then q p satisfies 

^ = q@q p + 9J, (6.2) 
where is the tetrad, linear in q and q p , defined in (4.11). Eliminating q p between (6.1) and (6.2) gives 

D 2 q Dq 



Dt 2 Dt 



©q+*P = q©q©q, (6.3) 
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which is not a completely closed because of the triplet of coefficients (1%, X%, s± in and the need to 
respect the Poisson equation. Through this, the vorticity is related to P and S by 



The associated eigenvectors of P and S project onto the the ortho- normal frame (£>, x, us x x) to yield 
the coefficients a, x, a p , c\ and c^. In this ortho-normal basis P has six components 



but for the present formulation, only a p , c\ and ci are required from P with a — G> ■ SCj and \ = 
(u> X x) ■ Su). 



Finally, recent developments in experimental and numerical capabilities also address the hybrid La- 
grangian and Eulerian descriptions of fluid dynamics while theoretical developments include the La- 
grangian averaged Navier-Stokes-alpha equations. The latter have been reviewed in Holm et al. (2005). 
The implications for Lagrangian averaging of the hybrid picture of rotating frames in ideal fluids pre- 
sented here will be discussed elsewhere but we note that the rotating frame representation may suggest 
a natural decomposition into fast and slow variables involving rapid rotations with slow modulations. 
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